Staggered orbital currents in the half-filled two-leg ladder 
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Using Abelian bosonization with a careful treatment of the Klein factors, we show that a certain 
phase of the half-filled two-leg ladder, previously identified as having spin-Peierls order, instead 
exhibits staggered orbital currents with no dimerization. 
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I. INTRODUCTION 

One of the most intriguing phases of strongly 
correlated electrons is kuawn variously as the "or- 
bital jaHti£Mcpmagnet,"BDa the "staggjered flux 
phase,"ljlaflM or the "d-density wave."E3y It is 
characterized by circulating currents which produce 
local magnetic moments aligned in an antiferromagnetic 
(staggered) way. As a consequence, time-reversal 
symmetry as well as translational and rotational sym- 
metries are spontaneously btipken. Another phase, the 
"circulating current phase, "t3 is somewhat similar, but 
does not break translational symmetry. These phases 
have received considerable attention lately, due to their 
possible relevance to the ijiseudogap region in the phase 
diagram ofpthe cuprates.liinl3 A recent neutron scattering 
experimentl; on underdoped YBa2Cu306.6 bas been 
interpretedliJ as evidence for these staggered orbital 
currents. 
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FIG. 1: Current fiow in the staggered fiux phase of the 
half-filled two-leg ladder. Reversing the currents gives the 
time-reversed state. 



In this paper we focus on the half-filled two-leg ladder, 
which is the simplest system that can support the stag- 
gered flux (SF) phase (see Fig. |l|). As the order param- 
eter of this phase breaks a discrete {Z2) symmetry, the 
possibility of true long-range order (LRO) of the currents 
is not a priori excluded at zero temperature in this one- 
dimensional system, in contrast to the situation for or- 
der parameters that break continuous symmetries, which 
causes their correlations to exhibit at most quasi-LRO 
with power-law decay. 

For weak interactions, the ladder can be treated using 
bosonization and the perturbative renormalization group 
(RG). For this case, the SF phase has been founcLin the 
phase diagram for spinless electrons at half-fillingJi3 Fur- 
thermore, away from half-filling, regions with dominant 
tendencies toward SFjjsi-dcring have been found both for 
spinlessEJ and spinfultZl electrons. Note that for general 
(i.e., incommensurate) fillings, true LRO of the currents 



is not possible due to the absence of Umklapp interac- 
tions (see Sec. ^ for a more detailed discussion). The 
results for the doped ladder with and without spin were 
summarized in Ref. [T8| which also investigated the effects 
of disorder. 

Here we are concerned with the SF phase for spinful 
electrons in a weakly interacting half-filled two-leg lad- 
der. In contrast to the other weak-interaction studies 
mentioned so far, in the approach used here the nearest- 
neighbor hopping parameters t_L and t along the rungs 
and legs, respectively, can be of the same order. We re- 
analyze the nature of a specific phase found in Ref. |l^, 
and demonstrate that this phase, previously identified to 
be of spin-Peierls (SP) type, actually exhibits staggered 
orbital currents with no dimerization, and therefore in 
fact is the SF phase. In contrast to the infinite-U half- 
filled case, where the constraint of no double occupancies 
makes the currents an unobservable gauge artifact, the 
LRO currents found here are observable. Furthermore, 
since all excitations are gapped, the SF phase at half- 
filling is insulating. 

We perform our calculations using Abelian bosoniza- 
tion, paying .careful attention to the Klein factors in this 
formalism.L^til As a check of our treatment, we also re- 
produce the identification of the CDW phase found in 
Ref. |l^. Furthermore, we show that our results, are con- 
sistent with those found for the doped ladder .EJ 

The paper is organized as follows. In Sees. || and [II 
we discuss the ladder model and its continuum limit and 
bosonized form, closely following the approach of Ref. |lj. 
In Sec. [V we define various local order parameters, de- 
rive their bosonized expressions, and calculate their ex- 
pectation values in the SF phase. The results are dis- 
cussed further in Sec. 0. Some of the technical details 
have been placed in two appendices. 



II. THE HALF-FILLED TWO-LEG LADDER 
AND ITS CONTINUUM LIMIT 



A. Kinetic energy 

We consider a two-leg ladder where the electrons can 
hop only between nearest-neighbor sites along the rungs 



and legs. The kinetic energy then reads 



^0 






l^{m+l)c£s{m)-t±y^ c\^{m)c2s{m)+R.c. 



(2.1) 

The operators c^sifn) and c^^{rn), respectively, annihi- 
late and create an electron on site ?7i = 1, . . . , A^ on leg 
£ = 1, 2 with spin s =t, i, and obey {cis{m), cl,^,{m')} = 
^u'Smm'^ss' , with all other anticommutators vanishing. 
Periodic (open) boundary conditions are used along (per- 
pendicular to) the leg direction. Introducing even and 
odd combinations 



1 



e/o,; 



V2 



{cis ± C2s 



(2.2) 



and Fourier transforming along the leg direction, the ki- 
netic energy becomes diagonal in momentum space, de- 
scribing two uncoupled bands with dispersion eo/o(fc) = 
— 2icosfca =F t±, where a is the lattice constant. Taking 
t_L positive, the even (odd) combination gives a bonding 
(antibonding) band. We consider a half-filled system and 
tj_ < 2<, in which case the Fermi level is at zero energy 
and crosses both bands, thus giving rise to four Fermi 
points zLkpc/o which satisfy kpc + kpo — tt/o. 

We will assume weak interactions and focus on the 
low-energy, long-wavelength properties of the model, so 
that we may linearize Hq around the Fermi points. It 
will be most convenient to work with a coordinate-space 
representation of the Hamiltonian. For this purpose we 
decompose the band operator cxs{m) (A = e,o) into a 
sum of left- and right-moving slowly varying (on the scale 
of the lattice constant) continuum fields, 



cxsim) 



'a e 



-ikp\x ^ 



V'LA.(a;)+e^'=-^^7/;flAs(a;)], (2.3) 



where x = ma. The linearized kinetic energy can then 
be written H^ = J dxTio i where 



Tio ^ -iVF ^ UpRXs^xIpRXs - i^lxsdx^LXs 
Xs 



(2.4) 



In this expression and throughout the paper it is under- 
stood that products of fermionic (and bosonic) opera- 
tors that may be evaluated at the same point are to be 
normal-ordered. The bare Fermi velocity vp is the same 
for both bands and is given by vp = a^J{2t)'^ — t'^. 



B. Interactions 

The continuum description of general, but weak, finite- 
ranged, spin-independent interactions, to leading order 
in the interaction strengths, was carefully discussed in 
Refs. 09, p2l and 03. One can restrict attention to terms 
which are both marginal (i.e., consisting of four-fermion 
interactions with no spatial derivatives) and nonchiral 
(i.e., containing two right-moving and two left-moving 
fermions). These terms can be classified according to 



whether they conserve momentum or not. The Hamilto- 
nian density for momentum-conserving terms reads 



LA/i 



Xfi 
+ fxfiJRXxJLfifi- fxtiJRXX- JLfitiJ- (2.5) 

Here / and b refer to forward and backward scattering 
processes, respectively, and 

S 

ss' 

where a^ , a^ , and a^ are the Pauli matrices. The follow- 
ing general symmetries hold: b'^^ = b'^^ and /^^ = /og, 
where v = p,a. To avoid double-counting, we set 
f^y^ = 0. At half-filling the model also has particle-hole 
symmetry, which implies h'^^ = fejf,,, leaving six indepen- 
dent couplings of this type. 

Half-filling also allows for non-momentum-conserving 
(i.e., Umklapp) terms. The Hamiltonian density for these 
interactions reads 

^f -E{<M4v4Ap-<^4v-^iA,.+H.c.}, (2.8) 

A^ 

where e = o and 6 = e. Here we have defined 

IpXii = E V'pasCss'V'Pms'' (2-9) 



'PA^ 



1 

2^^ 



'^^ 'lpPXs{<^0-)ss'i'Pti.s' , 



(2.10) 



where e = —10^. We may take u'^^ = u'^^ since Ipx^i = 
Ip^x and Ipxfi — —Ipp,x- The latter result also implies 
JpAA — 0, so that we can take u'^^ — 0. In addition, 
particle- hole symmetry gives u^^ = wg^, leaving three in- 
dependent Umklapp couplings. Thus a total of nine inde- 
pendent coupling constants must be taken into account 
in this model of the half-filled two-leg ladder. 



III. BOSONIZATION 

In the Abelian bosonization formalismpd'EllEJ the 
fermionic field operators ippxs can be expressecLin terms 



of dual Hermitian bosonic fields 
1 



ij 



PXs 



Xs and 8xs asEI 
Kxsexp[i{P(f>xs + dxs)], (3.1) 



where e is a short-distance cutoff, and P = R/L = ±1. 
The bosonic fields satisfy the commutation relations 



[(pXs{x),(l>X's'{x')] 
[^Xs{x),ex'Ax')] 



[exsix),ex's'{x')]^0, (3.2a) 
nr6xx'6ss'e{x-x'), (3.2b) 



the latter result written for e ^- 0. Here Q{x) is the 
Heaviside function. The long-wavelength normal-ordered 
fermionic densities can be expressed in terms of the 
bosonic fields as V'p^^^V'-PAs = dx{(l)\s + P9\s)/2tt. 

The Klein factors k\s commute with the bosonic fields, 
and satisfy 



{kXs,Kx's'} = 2SxX'Ss 



(3.3) 



Note that the Klein factors used here are Hermitian (in- 
stead of unitary), since we follow the common procedure 
of neglecting the number-changing jppperty of the Klein 
factors in the thermodynamic limit .£3 

Charge and spin operators are now defined as 

1 



(3.4a) 
(3.4b) 
with similar definitions of the 6 operators. We also define 






v2 



(3.5) 



where r — ± and v — p,a. Again, similar de finit ions 
apply to the 9 operators. Both Eqs. (3.4) and (B.J;) are 
unitary transformations, which implies that the commu- 
tation relations for the new sets of operators are of the 



same type as those in Eq. (3.2) 



Next we consider the bosonized form of the Hamilto- 



/(I) 



(2) 



nian density H — Hq + Ti-j + Ti-j , which is most suc- 
cinctly expressed in terms of the variables 
The kinetic-energy density reads 



and 



VF v^ 



^0 = ^ V [(d^A.)' + id.0r.)'] . (3.6) 

The momentum-conserving part of the interactions can 



be written Ti,} 



U\ 



(la) 



Til , where 



^/'"^ = i E ^- [(5x<^™)' - [d^Oruf] . (3.7) 
■/,^J with hp ^ 1, h„ ^ -1/4. 



Here Aru = K [K, 
Furthermore, 



T^f ' = - 



1 



(27re)2 



2r6^„ cos 29^ p cos 2<?!)+<, 



cos2(/.+<,(2&;:^cos2<?!>_^ + 2f /^';cos26'_^) 
cos 20_p(f6+ cos 2(/i_^ + h^^ cos 29 ^^ 



(3.8) 



with 6^ = &go ± 46§o and f = Kc-f'«oX'«oT'«oX- Finally, the 
bosonized form of the Umklapp interaction density reads 



Wi 



(2) 



■ cos 204 



(2716)2™"-^+" 

+2ul^ cos 20+^ -I- u+ cos 2(j)_a- 



SruPcos29^, 



FuZ, cos 29^ 



(3.9) 



with u± = <„ ± 4<^. 

Since the Hermitian operator F obeys F^ = / p \^g^ 
eigenvalues F = ±1. Furthermore, since [i?, F] = 0, H 
and F can be simultaneously diagonalized. 



IV. THE STAGGERED FLUX PHASE 



A. Pinned fields 



pf the one-loop RG equations for 
c3 that some of the couplings re- 



Numerical intcgratic 
the couplings showalj'c 
main small, while the others grow (sometimes after a sign 
change) and eventually diverge. The weak-coupling RG 
flow must be cut off before it leaves the regime of its per- 
turbative validity. The ratios of the diverging couplings 
at the cutoff scale are found to approach fixed constants 
in the limit of asymptotically small bare couplings, with 
different sets of ratios corresponding to different phases 
of the ladder. In the SF phase b^,^ and b'^^jS^e negligible, 
while the diverging couplings are given byllB 



fp 

J CO 



-r = -bP = -b" = 

A J CO ^CO A CO 



il, = -2< = 2< = <? > 0. 



(4.1) 



The resulting low-energy effective Hamiltonian can be 
mapped onto an S0(8) Gross-Neveu model, whose inte- 
grability may be exploited to extract the exact energies, 
degeneracies axid quantum numbers of all the low-energy 
excited states.E^ However, for our discussion, a semiclas- 
sical reasoning will suffice. Since the single coupling con- 
stant g flows toward large values, in the semiclassical 
ground state the bosonic fields in the Hamiltonian will 
be pinned to values which minimize the cosine interac- 
tion Ti) ' + Tij ■ Note that this argument would not 
be valid if the cosine interactions were to contain both 
the dual fields (j)-a and 9 -a, since then the uncertainty 
principle would forbid both fields to be pinned. However, 
0_CT disappears from the cosine interaction because b"^ is 
negligible and 6+ — u+ = 0. The pinned fields are then 
4>+p, 4'+cr, 9-p, and 9-a-- Since all four bosonic modes rv 
are pinned, the SF phase has no gapless excitations at 
half-filling. 

The possible solutions for the pinned fields are found 
by minimizing (F|_ff |F), where |F) is the eigenstate of F 
with eigenvalue F (these solutions will depend on F, but 
the p hysic s will of course not, as will be seen explicitly in 
Sec. |IVB|). There are infinitely many solutions for the 



pinned fields that minimize (F|iJ|F). However, this mul- 
titude of solutions is only apparent; taking into account 
the fact that the bosonic fields are not gauge- invariant, it 
can be shown thai there are only two physically distinct 
ground states.E^I The pinned-field configurations that we 
will use to specify these ground states are given in Table 



Ground state 


F 


(0+p) 


i^+a) 


{e-p) 


{0-<r) 


SF 1 
SF2 






■K 














SF 1 
SF 2 


— 1 


TV 








7r/2 
n/2 


7r/2 
7r/2 


CDW 1 
CDW 2 






■K 






n/2 
7r/2 






CDW 1 
CDW 2 


— i 




TV 










7r/2 
Tv/2 



TABLE I: The F-dependent pinned-field configurations used 
for the ground states in the SF and CDW phases (we choose 
(FlKeT^i^oTir) — i; see Appendix N). The two configurations 
listed here for a given ground state are physically equivalent, 
as can be seen from Table |l|. 



B. Order parameters 

In this subsection we explicitly show t hat the phase 
characterized by the couplings in Eq. (4.1) is not of 
spin-Poierls type with a (tt, tt) modulation in the kinetic 
energy,E2l but instead is the SF phase. We first define 
the relevant order parameters. The fundamental defini- 
tion of the current operator comes from interpreting the 
Heisenberg equation of motion for the number operator 



ne{m) = y^ cj^(m)cfa(TO) 



(4.2) 



as a discretized continuity equation. We will assume that 
the SF phase is a low-energy phase of a lattice Hamilto- 
nian whose interactions commute with ni{m). This is, 
e.g., the case for density-density and spin-exchange in- 
teractions. Then the components of the current operator 
take their conventional forms (see Fig. H) 



JUm 



) = it_La >^ C2s(m)cis(m) — H.i 



(4.3) 



jiim) = ita ^ c]!^(m + l)cfe(m) -H.c. . (4.4) 



Furthermore, the local kinetic-energy operator is 



kiini) = -i^ 4s ("^ + ^)cis(m) +11.' 



4 

(4.5) 



For completeness, in our discussion we also include the 
number operator n£(m) itself, since we will later check 
that our calculations reproduce the results for the CDW 
phase found in Ref. |lj. 

Next we outline the derivation of the bosonized expres- 
sions for these order parameters. It will be convenient to 
define an auxiliary operator, 

GeiiTijU^v) = y^\clg{m + u)cis{m) + V -JI.c. . (4.6) 

s 

m-1 m m+1 




jj(m-l) 



jj(m) 



FIG. 2: Currents as defined in E qs. (14. 3|) and (^.4|). Current 
conservation is expressed by Eq. (jA.13 



Then J£ (to) = itaGi{m, 1,-1), ki{m) — —tGi{m, 1,1), 
and ne{m) = Ge{m, 0,1). The continuum version of 

Ge{m,u,v) will contain products of type : V'pas(^ + 
ua)^jpi\isi{x) :, i.e., with the argument of the field op- 
erators differing by a lattice constant when u = 1 (here 
we have temporarily included the normal-ordering sym- 
bol explicitly) . One can safely Taylor-expand within the 
normal-ordering symbol to obtain : ^pxs{^)'4'P'\'s'{x) : 

+ua : dxippxs{^)i'P'ys'{x) : (note that due to the 
normal-ordering, all order parameters will be measured 
with respect to their values in the noninteracting ground 
state). For now, we only keep the zeroth-order term in 
the Taylor expansion, and comment briefly on higher- 
order terms later. We find 



PXs 



The expectation value of the normal-ordered long- 
wavelength density ippxs''l'P\s is zero. Bosonizing 
^pxs^-P^'^ produces exponentials containing the fields 



y+p, 



'■'+<7, 



and (p-a- Bosonizing ipp-^^ipp^s produces 



exponentials containing the fields 9-p, 9 



and I 



Since 



and 



are dual to the pinned fields 9_ 



and 9-cr: respectively, they will fluctuate strongly due 



to the uncertainty principle, and the expectation value 
of exponentials of these fields will vanish. We are there- 
fore left with a term which contains products of type 
'^p\s'^-P\s- Bosonizing this produces exponentials con- 
taining the four pinned fields, so this term will have a 
nonzero expectation value. 

Next consider j±{m). Its continuum expression only 

contains products of type iI'pxs'4'p\b ^-i^d 'ippxs'^-PXs- 
Thus only the latter product will contribute to the ex- 
pectation value of this operator. Note that in order to 
calculate j±{m) no Taylor expansion is necessary, since 
both fermion operators are taken at the same value of m 
from the outset. 

Using kpc+kpo = 7r/a and 2t cos kpoa = t± to simplify 
expressions, we find 



it^a\-iy''(F^i{x))+c.c., 



-zt^aH-lY+^^{F_^ix))+c.c., 



(4.8) 
(4.9) 



0±M) 

0X"^)) 

(fc,(m)) - la^t^ - (U/2)2(-l)^+"(Fi(a;)) -|-c.c.,(4.10) 



ai-lY+"\F,{x))+^ 



Here we have defined the operator 



JRos 



P^Lo 



:t 



Res 



(4.11) 



(4.12) 



The expectation value o f F^ jx) is inde pendent of x. It 
then follows from Eqs. ( |4.8| ) and (O) that if currents 
exist, they will flow around the plaquettes in a staggered 
pattern as shown in Fig. |^. Current conservation is ex- 
pressed by 



{Mm)) = {Mm-l)) + (-lY{u{m)) 



(4.13) 



(see Fig. H). Since expression ( |4.9| ) for {je{m)) only con- 
tains the zeroth-order term in the Taylor series expa nsion 
of the field operators, we conclude from Eq. ( [4.13| ) that 
the higher-order terms do not contribute to the plaquette 
currents. 

The bosonized expression for Fp [x) can be written 



PA-) = ^Xd[ 



27re 



AA '^-''s'^As 



As 
s(f)+„ - dxxi 



exp [i (0+p + s(f)+„ - dxxO-p - sdxxO-a)]- (4.14) 
Here we have defined deo = ^ ^oc = 1, s =Ti= ilj and 



a{p) 




(4.15) 



Let \n; T) = |n(r)) (g) jF) be a simultaneous eigenstate of 
H and F. The eigenstate |n(F)) lives in the Hilbert space 
where the bo sonic operators act, while |F) was introduced 
in Sec. IV A . We now consider a particular ground state, 
denoted by |0;F), and calculate the expectation value 
of Fp{x) in this state. First we insert the completeness 
relation ([A^) between the rightmost Klein factor and the 



Ground state 


r 


Fi 


F-i 


SF 1 
SF2 


±1 
±1 






4i 
-4i 


CDW 1 
CDW 2 


±1 
±1 


4 
-4 







TABLE 11 : Th e quantity Fp = 2ne(Fp(x))/B as calculated 
from Eq. (4.17), for the ground states in Table |. The physical 
properties of these states are seen to be independent of the 
"gauge" F, as they should be. 



exponential in Eq. (4.14) 
introducing (f>^p ~ (p^p — 
expression 

(exp[i(0+p-|-s0+<^ 



and use Eq. (A£). Upon 
4>+n) etc., we encounter the 



dx\t 



sdxxt 



,)]). (4.16) 



By construction, the pinned tilde-fields have zero expec- 
tation values. We also define 9+p = 9+p etc., for the fields 
dual to the pinned fields. As the Hamiltonian is invariant 
under a sign change of any of these tilde-fields, and their 
commutation relations are invariant under a combined 
sign change of any two dual fields, this expectation value 
is independent of s and dxx, as these variables can only 
be ±1. A qualitative estimate for this expectation value 
is calculated in Appendix M. Denoting the expectation 
value by B, we obtain 

B 



(^p(-))-£:E<f(---A3) 



As 



exp [^{{cj,+p) + s{4>+a) - dx-x{0-p) - -sdxxiO-a))]. (4.17) 

Here we have suppressed the F dependence of the expec- 
tation values appearing after the summation sign. This 
expression can now be evaluated for the ground-state 
configurations for the SF phase in Table | by inserting 
values for the pinned fields and using Eqs. (p]^) and 
( |A3| ) . The results are listed in Table ||. In the SF phase 
{F^i(x)) is nonzero and imaginary, which implies that 
the currents are nonzero. Explicitly, we find 

{j^{m)) ^2i-lY{Mm)) =T^^-8t^a\-ir, (4.18) 

where the upper (lower) sign refers to ground state SF 1 
(SF 2). Furthermore, {Fi{x)) vanishes identically, so that 
there is no modulation in neither {ki{m)) nor {n£{m)). 
We have also shown that the ficst order contribution to 
{ki{m)) is zero in the SF phase.Es 

Finally, we note that the ground-state degeneracy can 
be broken in a formal way by adding to the Hamilto- 
nian a term proportional to the order parameter. Thus, 
for the SF phase, one can let H ^ H — ?7i±(l), where 
r] is an infinitesimal constant. Depending on whether 
?7 ^ 0, ground state SF 1 or SF 2 will have the lower 
energy. The small imaginary-valued symmetry break- 
ing term perturbs the purely real-valued Hamiltonian, 



selecting a particular ground state which is intrinsically 
complex- valued with large imaginary components in the 
many-body amplitude. 



V. DISCUSSION 

As an additional check of our calculations, we have 
also reproduced the results for the CDW phase found 
in Ref. Ik 
are opposite to the ones given in Eq. 



In this phase, the signs of b^ b^^, and u^^ 



(|4.lD . The same 
bosonic fields are pinned as in the SF phase, but their 
expectation values are different. The CDW phase also 
has a two-fold degenerate ground state; the pinned-field 
configurations we have used are listed in Table |. The 
rest of the calculation is identical to the one presented in 
Sec. IV, including the calculation of B in Appendix |^. 
Our results for (Fp(a;)) for the CDW phase are summa- 
rized in Table ||. We find that {Fi{x)) is nonzero and 
real, so that {ni{m)) is modulated. This phase has no 
currents, since (F_i(a;)) = 0. 

It is perhaps worth commenting more explicitly on how 
the Klein factors affect the calculation of the expecta- 
tion val ues o f the various order parameters considered 
in Sec. IV B . T wo as pects are important. First, the 
(kas'^as) i^ Eq. (4.17) contribute relative signs to the 



various terms in the (A, s) summation. These signs are 
crucial for determining whether {Fp{x)) is nonzero, or if 
it instead vanishes identically due to cancellations. Sec- 
ond, if {Fp{x)) is nonzero, the fact that {kxsK^^) is purely 
imaginary affects whether {Fp{x)) is real or imaginary, 
which in turn determines whether the expectation value 
of a given order parameter that depends on {Fp{x)) will 
be nonzero; see Eqs. (^)-(4.11). 

In this paper we haye used the so-called "field- 
theoretic" bosonizationBJ We have also performed the 
calculations r-asing the more rigorous "constructive" 
bosonizationEll (however, we still neglect the number- 
changing prope rty o f the Klein factors). In the latter 
approach, Eq. ( 3.2b ) is replaced by [(j)\s{x),9ysf{x')] — 
i{TT/2)S\x'Sss'Sg'n{x — x'). Consequently, the anticommu- 
tation between right- and left-moving fermions with the 
same band and spin indices must now be taken care of 
by the Klein factors, which therefore acquire an addi- 
tional R/L index. As a result, 12 different products Ti 
of four Klein factors appear in the Hamiltonian. One 
must identify all relations between the Fj, as these re- 
lations put, restrictions on the permissible sets of eigen- 
values Fi.E3 Thus the treatment of Klein factors is more 
complicated than in the field-theoretic approach, where a 
single operator F appears. However, the final results for 
the expectation|-sAalues of the order parameters are found 
to be the samer3 

Our results imply that the SF phase occurs in the phase 
diagram of a weakly interacting general|-S0(5) invariant 
model on the half-filled two-leg ladder.Ej However, the 
basin of attraction of the SF phase is not restricted to 
have S0(5) symmetry. In fact, for all bare couplings 



studied in Ref. |l9|, including attractive interactions that 
break SO (5) symmetry, it was found that the RG flow 
goes to the SO (5) subspace, where the SF phase is one of 
the "attracting directions." It would be very interesting 
to undertake a complete exploration of the parameter 
space, to see if the SF phase could possibly be reached 
from purely repulsive off-site density-density interactions, 
supplemented by various spin-exchange interactions. 

Next, we discuss the possibility of SF order away from 
half-filling. For generic incommensurate fillings, Umk- 
lapp interactions are absent. Thus the total charge mode 
0+p will be gapless (making the system metallic), so 
that (exp (icjj^p)) wiH-janish. The currents will then 
only show quasi-LROO'ta'LZl Strictly speaking, the sys- 
tem is then no longer in the SF phase, but shows a dom- 
inant tendency toward SF ordering. On the other hand, 
for commensurateJUlings, higher-order Umklapp interac- 
tions are present ,E3Ej so that if these interactions are not 
irrelevant ,e3 (j>^p may be pinned, making true LRO possi- 
ble. These conclusions are consistent with those obtained 
from symmetry arguments: In the absence of Umklapp 
interactions, the Hamiltonian is invariant under the con- 



tinuous symmetry 



i+p + c (i.e., the constant 



can take arbitrary values), and pinning of cjj^p is for- 
bidden by the Mermin- Wagner theorem. This theorem 
no longer applies when Umklapp interactions are present, 
since then the symmetry is reduced to a discrete one (i.e., 
c can only take particular values). 

Finally, we show that our results for the half-filled SF 
phase are consistent with the results obtained for the 
doped ladder for generic incommensurate fillings. In Ta- 
ble II of Ref. |l8| the values of the three pinned fields in 
the phase with dominant tendency toward SF ordering 
are taken to be {(j)+a} — 7r/2, {6~ p) — 0, (O^a) = 0. Us- 
ing the SF couphngs in Eq. ([4.1D, and taking F = — 1, 
one sees that these expectation values of the pinned fields 



r(lb) 



(2) 



indeed minimize Hj , and also Hj at half-filling when 
{4>+p) is taken to be an odd multiple of tt/2. 
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APPENDIX A: MATRIX ELEMENTS OF 
PRODUCTS OF KLEIN FACTORS 



APPENDIX B: EXPECTATION VALUES OF 
EXPONENTIALS OF PINNED FIELDS 



The Hermitian operator F = KcT'^cX'^oT'^oX enters into 

bosonic expressions for Hj and W^ [Eqs. (3.8)-(3.9)]. 
Its eigenvalues are F = ±1, and the associated eigenstates 
|F) obey (F|F') = Jrr'- The completeness relation in the 
space spanned by these eigenstates is 



r=±i 



|F)(F|=/. 



(Al) 



We want to calculate various matrix elements of KcsHos, 
which appear in the expectation values of the order pa- 
rameters considered in Sec. IV B. We have 



Using f |F) = F|F), one obtains 

(fIkctKotIt) = r(F|KcxKox|r), 

(-FkoT'toTlr) = r(-F|KoXKox|F). 



(A2) 



(A3) 
(A4) 



The complex conjugate of Eq. (A4) can be rewritten as 
(F|kcT'«otI - r) = r(F|KcX'«oxl - r). Le tting F -^ -F in 
this relation, and comparing with Eq. (A4), shows that 
the off-diagonal matrix elements are zero; 



(-F|KcsKos|r) = 0. 



(A5) 



Next, consider the equation (F|f |F) — F. Anticommut- 
ing the two inner Klein factors and inserting Eq. (Al) 
gives 



J2 (r|KeTA^oT|r')(r'|KcxKox|r) = -F. 



(A6) 



=±1 



Using Eqs. (A5) and (AS), we obtain (F|KcsKos|r)^ 
— 1, i.e. 



(F|KosKos|r) = ±i. 



(A7) 



This result is consistent with k^sKos being an anti- 
Hermitian operator, thus having a purely imaginary ex- 
pectation value. 

In order to determine the matrix elements, one can, 
e.g., fix the two diagonal matrix elements for one of the 
spin directions. The other matri x el ements are then de- 
termined from Eqs. (A3) and (A5). In this paper we 
choose to set (F|KoxKot|F) = +i. 

Finally, we consider the space spanned by the states 
|n;F) defined in Sec. IV B. The completeness relation in 
this space reads 



J2 E kr)(n;F|=/. 
n r=±i 

We will also need the matrix element 



(A8) 



{0;T\Kxsn-^Jn;T') = (F|KAsACxJF')(0(F)|n(F')) 

= (F|KAs«Ajr)(5rr'<5„o. (A9) 



In this appendix, we calculate a qualitative estimate of 
the expectation value in Eq. (4.16) by employing a sim- 
plified treatment of the Hamiltonian, in which the pinned 
fields are expanded to quadratic order around their ex- 
pectation values. This gives (F|iJ|F) — J^riy ^rv with 
[Hr„, Hr'u'] = 0, where Hru is of Klein-Gordon form. 



1 

2^ 
1 

2^ 



Ti.-u — 

where vj^ = 24(7/ (tte'^) and 



uK{dJ+,f + -{dA+uf + w'0^, , (Bl) 






, (B2) 



u^j,l-w,r. i<-{^,^ m 



Furthermore, we have defined (jjj^p = (p+p — {4>+p) etc. 
for the pinned fields, and 9+p = 9+p etc. for the fields 
dual to them. Written in terms of these fields, Hr„ is 
independent of F (and also of whether we consider the 
SF or CDW phase). Clearly, H-i, can be obtained from 
H+i, by letting 6^,^ -^ 0_j^ and ip^i, -^ 9-v In addition. 



[dj+,{x),4^+,{x')] = [dj-,,ix),e-,,ix')]. (B4) 

Thus all Hamiltonians Hr^, are equivalent. It therefore 
suffices to consider, e.g., H+,y. We expand the fields as 



'2uJn 






dx9+^{x,t) 



'uK\ L 






[a^ 



9#0 



i(qx-ujgt) _ t p-i[qx- 



x-u,t)\ (gg) 



where usq = uj^q, and a+^q and a+y„ are canonical boson 

operators satisfying [a^^q.a'^^ ,] — 5qq/. These expan- 
sions give the correct equal-time commutation relations 
and equations of motion. The Hamiltonian can then be 
written on diagonal form. 



9#0 



(B7) 



with u!q — \Ju^q^ +uKw'^. 

Next we consider the ground-state expectation value 
of exp [ic0+j/(a;)], where c is an arbitrary c- number. Let 



$. 



<!> 



+1" 



(B8) 



where $+,y {^\.^) contains the annihilation (creation) 
part of (t>+u- The ground state expectation value can 
be written 



where 



exp 



'-[^+,.ix),<^l,{x)]), (B9) 



[^+,ix),^l^{x)]^^uK 



\~^ e 



-tq 



Wo 



(BIO) 



In actuality the coupling constant g is not constant up to 
arbitrarily high momenta. Rather, g is really a function 
(7(g), with g{q ^ 0) = g, and g{q -^ 00) = 0. It follows 
that u, K and w also become momentum-dependent, and 
ujq acquires an additional momentum dependence. Thus, 
a more correct expression for the commutator is 



[$+,(x),<I>^,(:r)] = ^^H(g)A-(g) 



-eq 



q>0 



UJq{q) 



(Bll) 



For simplicity, we will assume that there is a character- 
istic momentum cutoff 1/A such that for q <C 1/A, g{q) 
is well approximated by g, and for g ::!> 1/A, g{q) w 0. 
Multiplying the integrand by [e^^"* -I- (1 — e^'^')], we then 
approximate g{q) « g in the term containing e~ '^, and 
g{q) « in the term containing (1 — e^'^'). With A » e, 
this gives 



[$+.(x),$;,(x)]«^5]| 

g>0 



uK- 



-Aq 



p-e9 



-Ag 



l^[Ho{z)--Yo{z)] + lHA/e\ 



(B12) 



where z = Kw^jK/u. Here Hq{z) is a Struve function, 
and ^0(2) is a Bessel function of the second kind. 



Due to the equivalence of the Hamiltonians H^ 
have 



we 



o«C0 + p\ 



o'C0+a 



ice- 



JcO- 



(B13) 



It then follow s from Eq. ( p39D that the expectation value 
in Eq. (4.16) is independent of s and d^^^, as these vari- 
ables are restricted to be ±1. Denoting this expectation 
value by B, we find 



B^jexp{-^[Ho{z)-Yo{z)]}. (B14) 



The e in the prefactor cancels the 1/e in the prefactor 
of Eq. ( [4.17| ). In the unpinned hmit {g ^ 0), B fa 
^J 2Ag / TTv F — > 0, while in the limit of maximum pinning 
(.g -> .«^), i? « i exp (-iy^iS) ^ ^. 
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